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Abstract 

The gauge invariance of the muonium-antimuonium (MM) oscillation time scale is explic- 
itly demonstrated in the Standard Model modified only by the inclusion of singlet right-handed 
neutrinos and allowing for general renormalizable interactions. The see-saw mechanism is ex- 
ploited resulting in three light Majorana neutrinos and three heavy Majorana neutrinos with 
mass scale Mr » Mw- The leading order matrix element contribution to the MM oscilla- 
tion process is computed in gauge and shown to be ^ independent thereby establishing the 
gauge invariance to this order. Present experimental limits resulting from the non-observation 
of the oscillation process sets a lower limit on Mr roughly of order 600 GeV. 

1 Introduction 

Muonium (M) is the Coulombic bound state of an electron and an antimuon {e~ ji^), while antimuo- 
nium (M) is the Coulombic bound state of a positron and a muon {e^jT). It was suggested roughly 
50 years agoIH that there may be a spontaneous conversion between muonium and antimuonium 
which would violate the individul electron and muon number conservation laws by two units. 
Such a muonium-antimuonium oscillation is totally forbidden within the Standard Model. Hence, 
its observation will be a clear signal of physics beyond the Standard Model. Since the initial 
suggestion, experimental searches have been conductedEl-lH and a variety of theoretical models 
have been proposed which could give rise to such a muonium-antimuonium conversion. These 
include interactions which could be mediated by (a) a doubly charged Higgs boson A"*""*"© |7l, 
which is contained in a left-right symmetric model, (b) massive Majorana neutrinos ISIHI, or (c) 
the T-sneutrino in an R-parity violation supersymmetric model ifTOl . 
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In this paper, we focus on a modified Standard Model whicli includes singlet right-handed neu- 
trinos. There is now compelling evidence of the existence of neutrino oscillations from the exper- 
imental study of atmospheric and solar neutrinos ifTTI - ifTSl . That implies nonzero neutrino masses 
and mixing matrix elements. The size and nature of the neutrino mass and the associated mixing 
is still an open question subject to experimental determination and theoretical speculation fl6\ - 
ifTSl . One simple neutrino mass model is obtained by modifying the Standard Model by including 
singlet right-handed neutrinos and allowing for a general mass matrix for neutrinos. Left-handed 
neutrinos along with their charged leptonic partners are components of SU {2)l doublets and ex- 
perience the weak interaction while any right-handed neutrinos are completely neutral under the 
Standard Model gauge group. The see-saw mechanism |[T9l - |[23l provides a natural explanation 
of the smallness of the three hght Majorana neutrino masses, while ensuring that the other three 
Majorana neutrinos are heavy. Such a model could also lead to the muonium-antimuonium oscil- 
lation process. In order for there to be a nontrivial mixing between muonium and antimuonium, 
the individual electron and muon number conservation must be violated. Such a situation results 
provided the neutrinos are massive particles which mix amongst the various generations. This 
criterion can be met by the modified Standard Model and the e~fj.'^ and e^/u' states could indeed 
mix. 

2 Neutrino masses and mixings 

The leptonic Yukawa interactions with the Higgs scalar doublet in the modified Standard Model 
take the form 

Here and vfj are respectively the charged lepton and its associated neutrino partner of the 
5 U(2)l doublet, while vj^'j is the right-handed neutrino singlet. The superscript zero indicates 
weak interaction eigenstates so that the leptonic charged current interaction is 

^".-4^"^'ZO.vS + //.C. (2) 

After spontaneous symmetry breaking, the mass term for the charged leptons takes the form 

3 

J-mass -- 2_i ^^Ra^ab'^Lb + ^La^ba'^Rb^ 
a,b=\ 



2 



where is a 3 x 3 mass matrix. To diagonalize this matrix, one peri'orms the biunitary transfor- 
mation 

'"'=^M,a,(^'')^ (4) 

where and are 3 x 3 unitary matrices and m^^.^^ is a diagonal matrix whose entries are the 
charged lepton masses. To implement this basis change, the charged lepton fields participating in 
the weak interaction are rewritten in terms of the mass diagonal fields as 



b=l 



b=l 



So doing the mass term reads 



a=l 



A general neutrino mass term resulting from renormaUzable interactions takes the form 



(6) 



£V ^ _!/ . (0).^ (0) \ 



{m^Y 
rrp 



( (0) ^ 



+ H.C. 



(7) 



Note that the upper left 3x3 block in the neutrino mass matrix is set to zero. This block matrix 
involves only left-handed neutrinos and in the (modified) Standard Model its generation requires 
a nonrenormaUzable mass dimension-five operator. Consequently such a term will be ignored. 

For three generations of neutrinos, the six mass eigenvalues, ruyA, are obtained from the diag- 
onalization of the 6x6 matrix 

{ (m^)^ 1 



(8) 



Since Af is symmetric, it can be diagonalized by a single unitary 6x6 matrix, JJ, as 



(9) 



This diagonalization is implemented via the basis change on the original neutrino fields organized 
as the 6 dimensional column vector 



( (0) ^ 



^^R 



\ ^R ) 



(10) 



to the new neutrino fields defined as 



(11) 
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where 



Vl 


, Nr = 













(12) 



The neutrino mass term then takes the form 

^6 6 

■Cmass = ~2 ^ myA\v\CvA + vJCv^] ^ - ^ niyj^VA-, (13) 
A=\ A=l 

where niyA are the Majorana neutrino masses. 

Since a nonzero Majorana mass matrix does not require S U{2)l x U(l) symmetry break- 
ing, it is naturally characterized by a much larger scale, Mr, than the elements of the matrix 
whose nontrivial values do require SU{2)i xU{\) symmetry breaking and are thus expected to be 
somewhere of the order of the charged lepton mass to the W mass. Thus one can take the elements 
of nP , characterized by a scale wid, to be much less than Mr, the scale of the elements of m^. One 
then finds on diagonalization of the 6x6 neutrino mass matrix that three of the eigenvalues are 
crudely given by 

2 

in 

'"I'a ~ ITT ^ ^D, a = 1,2, 3, (14) 

Mr 

while the other three eigenvalues are roughly 

rUyi-MR, /-4,5,6. (15) 

This constitutes the so called see-saw mechanism |fT9l - 1|23]| and provides a natural explanation of 
the smallness of the three hght neutrino masses. Moreover, the elements of the mixing matrix are 
characterized by an Mr mass dependence 

Uah-0{1), a,b= 1,2,3 
Uij-^OiD, i,j = 4,5,6 

Uia ~ Uai ~ (?(^), a = 1,2,3, i = 4,5,6. (16) 
Mr 

Since the charged lepton mixing matrix is independent of Mr, one finds that elements of the mixing 
matrix appearing in the charged current have the Mr mass dependence 

Vah-0{\), a,b =\,2,3 

Vai~Oi^), a -1,2,3, / = 4,5,6 (17) 
Mr 

Inserting the transformations (5) and (11) in the interaction terms (1) and (2), and taking 
into account the mass matrix transformations (4) and (9), one obtains the explicit interactions of 
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charged bosons with the leptons in their mass diagonal basis as 

3 6 3 6 

= -^W->^ 2] 2] fLay^VaAVA ' ^W^^ J] Z ^aKaT/l. (18) 
a=i A=l v2 a=l A=l 

Xf„j = (p" y y (mia ^ - m^A^-^) VA 

-^^-^^ Z Z ^^^«A K-^^ - 4 (19) 

where 

3 

V«A = Y}-Al\bUhA (20) 
Note that the mixing matrix Vca satisfies the identities 1231: 

6 

Yj ^^aVIa = 5ab (21) 
A=l 

6 

2] ^mVm'^^va = (22) 

A=l 

Identity (21) stems from the unitarity of matrices Ai and U , while identity (22) is a consequence 
of the particular form of the neutrino mass matrix. In particular, it requires the Majorana mass 
term of the left-handed neutrinos be set to zero. A detailed proof of this later identity is provided 
in the Appendix. 

3 The gauge invariant T-matrix elements 

The lowest order Feynman diagrams accounting for muonium and antimuonium mixing are dis- 
played in Fig. 1. We shall consistently employ the gauge. The gauge invariance of the T-matrix 
element will be demonstrated by establishing its ^ independence. In Fig. 1, there are two neutrinos 
in the intermediate state for each graph while every wavy line represents either a W boson or an 
gauge charged erstwhile Nambu-Goldstone boson. 
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(d) 



Fig. 1 Feynman graphs contributing to the muonium-antimuonium 
mixing. Each wavy line is either a W boson or an gauge charged Nambu-Goldstone boson. 



Note that in unitary gauge ( ^ ^ oo), the W boson propagator takes the form 2 — -[gfiv - 
^^]. A theory with such a propagator has very bad power counting convergence properties. As it 
turns out, the unitary gauge power counting divergent pieces in the W vector box diagrams vanish, 
as they must, after application of properties (21) and (22). Hence, when we calculate the T-matrix 
elements in gauge, we will also apply properties (21) and (22) to establish the cancellation of 
the various terms in this case. 

As it turns out, graph (a) gives the same contribution as (b), as do graphs (c) and (d). Hence, 
we need only discuss the gauge invariant T-matrix elements of the graphs (a) and (c). Fig. 2 details 
expUcitly the 4 separate graphs which are represented by the single graph in Fig. 1. 
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J" 




(al) 

y" e 





(a2) 



(a4) 



Fig.2 Feynman graphs of type (a) in gauge. 



A straightforward application of the gauge Feynman rules 11251 to the above graphs yields 
the T-matrix elements 



6 6 



■f 

Jo 



dt 



XaXb 



+1)2 I t 



A=l B=\ 

2(^-1) ^ ^ (^-1)2 



(? + xaX? + xbX? + 1)2 I" f + ^ 4(f + ^)2 

6 6 



(23) 



7"a2 - Ta3 



1 ^75 
2 



Jo 



^- 1 2 

f + -2 r 



7"a4 - 



647r2M2, 



a + Xa)(? + XB)(f + l)(f + f) ^ 4(f + ^) 
1 - 1 - ^ ^ 



(24) 



A=l B=l 



r 



dt 



(25) 



{t + XA){t + XB){t + ^)'^ 4 

where /i(3) = il(pj,sj) , /i(4) = il(p4,S4) , e(l) - e(p\,si) and e(2) - e(p2,S2) are the spinors 
of the muons and electrons and xa - , A - I, 6. Note that in obtaining these results, we 
akeady applied properties (21) and (22) to eliminate various self-cancelling terms. As such the 
integrals in (25)-(27) are finite even in the ^ ^ oo limit. 

In order to discuss the ^ dependence in a manifest way, we rewrite these T-matrix elements as 

T - matrix elements in gauge 



6 6 -j^ 

= Z Z dtnxA, XB, t) ■ (hit) ■ - — - igit) ■ — + m) 

Ta2 = YjYu I dtJl(XA,XBj)-{-h(t)-——-^+g(t)-—) 

^«3 = yy dt^(xA,XB,t)-{-h{t)- zr^ + g(t) ■ —-) 

6 6 poo 

TaA = YjYj I dtJl{XA,XB,t)-(h{t)-——^) 




where 

^{xA,XB,t) = - ^J,^, [iU(3)y;.^^g(2)][^(4)y'^2:^g(l)](V^V;^)(V^BV;g) 

(26) 

(/ + XA){t + Xb) 

with 

/J(0=4 , ^^^^ = 771 '^"^ /(0= ^^^^^2^ (27) 
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Note that the ^^^^^2 terms from the second and third graphs totally cancel against the ones from 
the first and fourth graphs ,while the terms from the second and third graphs exactly cancel the 
one from the first graph. All ^ dependent contributions thus vanish and the only remaining piece 
is the term containing f(t) from the first graph, which is ^ independent. Hence, we have the gauge 
invariant T-matrix element for graphs of type (a) 



Ta = 



1 - 1 - 



ut/i lYi^ - - A=1S=1 

1 + 2f + j 
, dt T- 

it + XA)it + XB){t + 1)2 

GlMl 



Jo 



^ '^D«(3)r;.(l - 75)^(2)] [;i(4)/(l - 75)^(1)] 



87r2 

6 

+ Yj (VO(VKb)7^(^a,xb) 

A,B=\\A+B 

Here we have introduced the Fermi scale 



Y}y,Av:jfs{XA) 



A=\ 



V2 ml 

along with the Inami-Lim|26| function 



We have also defined 



(28) 



(29) 



x^-\\x^+Ax 3x3 
4(1 - xY 2(1 - xy 



T{Xa, Xb) = XAXBi"^^^-^ AfB)\ ^ j.^^^^ ^21) 

^ xa- Xb ' 

with 

4(1 - xY 4 (1 - x) 



In a similar manner, the graph (c) in Fig.l represents the following four graphs: 
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e e e e 



(c3) (c4) 
Fig.3 Feynman graphs of type(c) 

The T-matrix elements of the above four graphs are 

T - matrix elements in gauge 
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V 



6 6 



6 6 ^oo ^ 



where 



mA • triB xaXb 



{t + XA){t + Xb) 



(33) 



and 



+t + 2 



f{XA,XB,t)= (34) 

In a similar fashion to the case for the graphs of type (a), all the ^ dependent terms again 
cancel against each other leaving only the ^ independent f{xA,XB,t) term. Thus the type (c) T- 
matrix element is gauge invariant and is given by 

^"F^.r....^.. ................. ^2 V...* .2 



= ^V^bu(3)/(1 - r5)e(2)][ju(4)7^(l - 75)^1)] XCV)' TSKb) 



2 

If xa - Xb, the relevant integral is 



r J 4t^xAXBit^2) I 

Jo 1(? + xa)(? + xb)(?+1)2/ 



2 Jo (t + XA){t + XB){t + l)^ 



i-^A - 4)XA ^ (x^ - 3a;^ + 4a.i + 4)xa 
1 + ^7 ^3 ^^^^ 

{xA - \Y 2{xA - \y 



while for xa + xb, it takes the form 



K(xA, Xb) - 



^J XAXB r°°jj> 4t + XAXB{t + 2) 

2 Jo (t + XA)(t + XB)(t + 1)' 

L(xA, Xb) - L(xb, xa) 

yXAXB (37) 

Xa - Xb 
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with 



, , , 4 - xaXb xa(2xb - xaXb - 4) 

L(xa, xb) = — + In xa (38) 

2{xa - 1) 2{xa - ly 



The T-matrix element of graph (c) is thus secured as 



Tc - -^[;i(3)/(l - 75 M2)] 1^(4)7^(1 - rsMl)] 
6 6 



{Y}y,Av:Afi{XA)+ (V)'(Kb)'^(^a,xb)] 09) 



A=l A,B=1;A#B 



Combining the various contributions, the T-matrix element can be reproduced using the gauge 
invariant effective Lagrangian given by: 



G 



MM 



|^yj(l _ y5)e][jU7^(l - 75)e] 



(40) 



where 



^MM 

V2 



6 



2](VKa)''^(^a)+ 2] (VKa)(VKb)7^(^a,^b) 

A=l A,B=l;A5tB 

6 1 

2]( VKa)''^(^a) - 2] (V)'(Kb)'^Ua,X5) 
A,B=l;A7tB 

6 

2](VKa)'('^(^a)-/(xa)) 



16;r2 

6 



A=l 



+ 2 (( VKa)( VKb)7^(^a, xb) - ( V)'(^:b)'^(^a, xb)) 



A,B=\\A+B 



(41) 



4 Limit on M 



Muonium (antimuonium) is a nonrelativistic Coulombic bound state of an electron and an anti- 
muon (positron and muon). The nontrivial mixing between the muonium ( |M > ) and antimuo- 
nium (|M >) states is encapsulated in the effective Lagrangian of Eq. (001) and leads to the mass 
diagonal states given by the linear combinations 



1 



\{\ + e)\M > ±{\ - e)\M >] 



where 



V2(l + leP) 

^M-MM - ^M-MM 

f- — - 



(42) 



(43) 
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< Ml - r d^r£,ff\M > <M\- f d^r£,ff\M > 

Mmm- , i- , Mmm - , i- - (44) 

V< M\M >< M\M > V< >< M\M > 

Since the neutrino sector is expected to be CP violating, these will be independent, complex 
matrix elements. If the neutrino sector conserves CP, with \M > and \M > CP conjugate states, 
then Mf^f!^ = Mi;^}^^ and e = 0. In general, the magnitude of the mass splitting between the two 
mass eigenstates is 



|AM| = 2 



Re yjMMMMMM 



(45) 



Since muonium and antimuonium are linear combinations of the mass diagonal states, an initially 
prepared muonium or antimuonium state will undergo oscillations into one another as a function 
of time. The muonium-antimuonium oscillation time scale, r^^, is given by 

^ \AM\. (46) 



'''mm 

We would like to evaluate \AM\ in the nonrelativistic limit. A nonrelativistic reduction of the 
effective Lagrangian of Eq. (l40l) produces the local, complex effective potential 

y.//(r) ^ S^^S\r) (47) 

Taking the muonium (anitmuonium) to be in their respective Coulombic ground states, 0ioo(r) - 
g-riaMM where a^u = — ■ — is the muonium Bohr radius with ntrpd - — r-^ - the re- 

duced mass of muonium, it follows that 



1 



: J d'rcl>\Qoir)\ReV,ff(rMr)m 



'MM 

16 — — |^ioo(U)l p n — (48) 

MM 



Thus we secure an oscillation time scale 



1 ^ 16|/?eG^Ml^3^3 (49) 



Tmm ^ V2 

The present experimental IrmitlH on the non-observation of muonium-antimuonium oscilla- 
tion translates into the bound \ReGf^,^\ < 3.0 x IQ'^Gp where Gp ^ 1.16 x \0'^GeV~^ is the 
Fermi scale. This limit can then be used to construct a crude lower bound on Mr. For the case 
when the neutrino masses arise from a see-saw mechanism and taking to be of order Mw, the 
Mr dependence of is obtained from Eq. (|4T]) as: 

GlMt Mr 

Casel: IReG^ul ^-^In— ^, A = 1, 2, 3, B - 1, 2, 3 

M^ Mw 
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Case!: l^eG^wl ^-^In— , A = 4, 5, 6, - 4, 5, 6 

GlM% Mr 

Case3: \ReGMM\ —r^ln—, A = 1, 2, 3, B - 4, 5, 6 (50) 

Case 1 and case 2 give the same order Mr dependence, while case 3 is suppressed by an 
additional factor of Hence, the term In ^ gives the dominant contribution. We then 
roughly calculate a bound of Mr as 



Gt-Mt, Mr -, 
F W i„ ^ T A w in-3 



2 fl^4 

, In— < 3.0x IO^^Gf (51) 
Ml 

which has also been obtained in reference 191. Using My^/ ~ 80.4 GeV and Gp - 1.166 x 
IQr^GeV'^, we finally secure 

Mr > 6 X IQ^GeV (52) 

Note that this is just a rough estimate since we are retaining only the dependence on Mr while 
neglecting all numerical dependence on the mixing angles and CP violating phases in VaA- 
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Appendix: Proof of Identity (22) 



Using the definition of the mixing matrix VaA = 'Lc=ii^L )ac^cA, one can write 

6 6 f 3 \ f 3 ^ 



A=l 



Y^VaAVbAtnyA = 2 Y^iAi^\cUcA ' ^u'^A^^^U^A 



A=l \c=l 
3 3 



f 6 



c=l «i=l VA=1 

where myA are the diagonal elements of matrix M'^.^^, 

myA = (M^diag^AA. 



(53) 



(54) 



Consequently, we can express H^^^ UcAmyAUdA as a product of matrices and equation (53) takes 
the form 



3 3 



2 yaAVbAmyA = 2 X^'^^')- (^^Lg^O,, (^41')^ 



(55) 



A=l c=l (i=l 

fv _ ttTmvi 



Using equation (9), M^.^^ = f/^ M^f/, it follows that 



Substituting this result back into eq. (55) then gives 

6 3 3 



(56) 



J] VaAVbAmyA - Tu^^~L)ac(Mncd(Al\d 



A=l 



c=l d=l 



where 



(m^)' 



(57) 



(58) 



Since c and d both run from 1 to 3, (M'^*)cd are the elements of the upper left 3x3 block of matrix 
(58), which is zero. Hence, we secure the identity 



^ VaAVbAmyA - 



(59) 



A=l 
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